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We systematically analyze the influence of the superconducting gap symmetry and the electronic 
structure on the dynamical spin susceptibility in superconducting Naa;Co02-J/H20 within a three 
diflerent models: the single aig-band model with nearest-neighbor hoppings, the realistic three-band 
t29-model with, and without e'g pockets present at the Fermi surface. We show that the magnetic 
response in the normal state is dominated by the incommensurate antiferromagnetic spin density 
wave fluctuations at large momenta in agreement with experimental temperature dependence of 
the spin-lattice relaxation rate. Also, we demonstrate that the presence or the absence of the e'g- 
pockets at the Fermi surface does not affect significantly this conclusion. In the superconducting 
state our results for d^2_y2- or dxy-wave symmetries of the superconducting order parameter are 
consistent with experimental data and exclude nodeless d^2_y2 + id^y-wave symmetry. We further 
point out that the spin-resonance peak proposed earlier is improbable for the realistic band structure 
of Naa;Co02-J/H2 0. Moreover, even if present the resonance peak is confined to the antiferromagnetic 
wave vector and disappears away from it. 

PACS numbers: 74.70.-b; 75.40.Gb; 74.20.Rp; 74.25. Jb 



I. INTRODUCTION 

The spin dynamics in unconventional non s-wave su- 
perconductors is of fundamental interest due to its inter- 
esting and peculiar properties. This includes a non-trivial 
behavior of the magnetic part of the Knight shift in the 
spin-triplet superconductors^, as well as an emergence of 
the so-called resonance peak observed in superconduct- 
ing layered cuprates^ which possesses spin-singlet (1^2 _y2- 
wave order parameter symmetry. Furthermore, magnetic 
excitations are also often considered as a possible glue 
for the Cooper-pairing in a number of hcavy-fcrmion and 
transition metal oxides compounds. 

An analysis of the feedback effect of superconduc- 
tivity on the magnetic spin susceptibility can be used 
to determine the symmetry of the superconducting or- 
der parameter. This is of particular significance for re- 
cently discovered water intercalated sodium cobaltate 
superconductor—, Naa:Co02 • yH20, where the origin of 
superconductivity as well as an underlying symmetry of 
the superconducting order parameter is currently under 
debate. The studies of the specific heat^i^>^ii and the /xSR 
measurements of a magnetic penetration depth^ have re- 
vealed a line of nodes in the superconducting gap func- 
tion Ale. Similar conclusion has been made based on the 
measurements of the spin-lattice relaxation rate 1/TiT 
by means of Nuclear Quadrupolc Resonance (NQR), 
where absence of the characteristic Hebel-Slichtcr peak 
and power-law decrease upon decreasing temperature has 
been observed^ i^°'^^i^^'^'^ . Simultaneously, the develop- 
ing of the strong antiferromagnetic (AFM) fluctuations 
above superconducting transition temperature, Tc, have 
been found. At the same time, early reports on the 



Knight shift's temperature dependence, K{T), have sug- 
gested a spin-triplet symmetry of the superconducting 
gap ^^i^^ . In these Nuclear Magnetic Resonance (NMR) 
experiments, K(T) was shown to be anisotropic for exter- 
nal magnetic field applied parallel or perpendicular to the 
a6-planc. In particular, Kc{T) component has not shown 
a substantial decrease below Tc- This behavior has been 
interpreted in favor of the odd-parity Cooper-pairing in 
sodium cobaltate a-'^^'-'^'^i-'^^i^^i^^'^"'^ . However, the most re- 
cent NMR experiments with higher precision have found 
a reduction of both Knight shift components as a func- 
tion of temperature for T < T,^^^. These experiments 
points towards spin-singlet Cooper-pairing. 

From the group-theoretical analysis the even- 
parity symmetries of the lowest harmonics 
for the triangular lattice are classified accord- 
ing to s-wave (Ak = Aq), extended-s-wave 
(Ak = 2/3Ao[cosfey + 2cos{k^Vi/2)cos{ky/2)]), 
dx2_y2-wave (Ak = Aq [cos fc.y — cos (fc,V3/2)cos(fcy/2)]), 
dxy-v/avc (Ak — Ao[V3 sin (fca;V^/2) sin (fcj^/2)]); and 
dx2_y2 -\- idxy-'^&YB representations^^. For both dx2_y2- 
wave and dxy-wa,Ye symmetries Ak has line of nodes 
at the Fermi surface. Moreover, the time-reversal 
symmetry is broken for dx2_y2 -\- id^y-'wave state. 

For the pure trigonal symmetry of the Co02-plane, 
all three d-wave states are degenerate. However, due 
to the absence of nodes dx2_y2 -\- idxy-'wave seems to 
be most energetically favorable. Until now, a break- 
ing of time-reversal symmetry has not been observed 
in experimen t^^i^^ . Generally, the combined influence 
of the impurities and some competing instabilities, such 
as Cooper-pairing in a secondary channel as well as the 
lattice symmetry breaking, can lift the degeneracy be- 
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twcen these three d-wave competing ground states^l. 
This may indeed be the case for sodium cobaltates 
where Na arrangement introduces disorder at a; = 0.33 
concentration^. More sophisticated theories, involving 
multi-orbital model for sodium cobaltates, suggest two 
different gap symmetries (one of which is dr^2_y2 -|- idxy) 
for two different Fermi surface topological^. 

Obviously, there is still a controversy on the symmetry 
of the superconducting order parameter in sodium cobal- 
tates. In present study we systematically analyze the in- 
fluence of the superconducting (SC) gap symmetry and 
the electronic structure on the dynamical spin suscep- 
tibility in Naa:Co02-?/H20. In particular, assuming spin 
singlet s-wave and d-wave symmetries of the supercon- 
ducting order parameter we have calculated the real and 
the imaginary part of the magnetic response as a function 
of the momentum, temperature and frequency. We de- 
duce the characteristic temperature dependencies of the 
Knight shift and spin-lattice relaxation rate. Further- 
more, we have studied the feedback of the superconduct- 
ing order parameter on the frequency dependence of the 
imaginary part of the spin susceptibility. We investigate 
the role played by the details of the electronic structure 
of Na2;Co02-yIl20 and, in particular, the changes of the 
Fermi surface (FS) topology induced by the multi-orbital 
effects. 

Structurally, a parent compound, Na2;Co02, has a 
quasi-two-dimensional structure with Co ions in the 
C0O2 layers forming a triangular lattice. Na ions re- 
side between these layers and donate x electrons to the 
partially filled Co-d{t2g) orbital. Apart from doping, Na 
ions also induce structural ordering at higher doping con- 
centrations [x > 0.5) where superconductivity does not 
occur. Due to the presence of a trigonal crystalline elec- 
tric field (CEF), the t2g level splits into the higher ly- 
ing aig singlet and the two lower lying e'g states. The 
ab-initio band structure calculations within a Local Den- 
sity Approximation (LDA) predict Naa;Co02 to have a 
large Fermi surface with mainly aig character and six 
hole pockets of mostly e'g character—. At the same time, 
surface sensitive Angle-Resolved Photo-Emission Spec- 
troscopy (ARPES)^!^!^ reveals a doping dependent 
evolution of the Fermi surface, which shows no sign of the 
e'g hole pockets for 0.3 < a; < 0.8. Instead, the observed 
Fermi surface is centered around the F point and has 
mostly aig character. It has been argued that such an ef- 
fect may arise due to strong electronic correlations^^j^^i^, 
however, no consensus in the literature has been reached 
yet (see e.gM:^^). 

In Na2;Co02-j/H20 due to the water intercalation the 
inter-layer C0O2 distance becomes larger and, thus, the 
material becomes more two-dimensional leading to a sub- 
stantial decrease of the bilayer splitting. However, little 
is known about the particular changes in the electronic 
structure and the energy splitting between aig and e'g 
levels. 

In order to take into account the multi-orbital effects 
wc analyze the effect of superconductivity for the three 



different cases: the single-band (aig) model with nearest- 
neighbor hoppings, the realistic three-band (^29) model 
with, and without six e'g pockets at the FS. 



II. 



aig- 



BAND MODEL 



We first consider the simple aig-band model, repre- 
sented by a two-dimensional Hubbard Hamiltonian on 
the triangular lattice: 



H 



k,o- 



flkcr 
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where ni„ — a\^ai„, aia {a\^) is the annihilation (cre- 
ation) operator for the aig hole at the Co site i with spin 
(J. Here, ek = 2t[cosky + 2 cos {kx^/i/ 2) con {ky/ 2)] — fi, 
t=0.123 eV is the nearest-neighbor hopping integral, and 
/i is the chemical potential which has been calculated 
sclf-consistcntly for x = 0.33. The energy dispersion, Sk, 
along the principal directions of the hexagonal Brillouin 
zone (BZ) and the corresponding Fermi surface are shown 
in Fig.[TJb) and in Fig. [21 respectively. Here, F = (0,0), 
K = (0,2/3), and M = (1/2^3,1/2) [in units of 27r/a] 
denote the symmetry points of the first BZ. Later, coor- 
dinates of the wave vectors will be given in units of 27r/a 
with a being the in-plane lattice constant. 

To calculate the dynamical spin susceptibility, we em- 
ploy the Random Phase Approximation (RPA) which 
gives 



Xo(q, iuJrn) 
1 - C^Xo(q, ii^m) ' 
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where xoi^.ii^m) is the BCS Lindhard susceptibility 
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with 



^^ e^e^+.,+Rc(A^A^^j being the BCS coher- 
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ence factors. Here, Um are the Matsubara frequencies, 
f{E) is the Fermi function, and i?k = + |Akp. 

In Fig. [ija) we show both the bare and the RPA mag- 
netic susceptibility in the normal state at w = 5 meV and 
U = 0.25 cV. One immediately notices that the magnetic 
response is dominated by the scattering at the incommen- 
surate wave vector, Qsdw = (0,0.598) ~ (0,3/5). The 
value of Im[x(q, w)] at the commensurate wave vector, 
Qafm = {(0, 2/3), (l/-\/3, 1/3)}, appears to be much 
smaller. There is also another incommensurate wave vec- 
tor present, Qsdw- presence of a set of incommen- 
surate wave vectors with substantial magnitude of mag- 
netic scattering shows a tendency of the itinerant elec- 
trons on the triangular lattice towards spin density wave 
(SDW) instability. 
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FIG. 1; (Color online) Calculated results for the aig-band 
model, (a) q-dependence of Im[xo(q, t^^)] and Im[xRPA{<l, i^)] 
at w = 5 meV in the normal (non-SC) phase. The scatter- 
ing wave vectors Qafm, Qsdw, and Qsuw a-re denoted by 
the arrows, (b) Calculated aig-band dispersion, where the 
horizontal (green) line stands for the chemical potential. The 
panels (c)-(e) show imaginary and real parts of xo, and imag- 
inary part of xrpa at q = Qafm in non-SC phase and in SC 
phase with various superconducting order parameter symme- 
tries. The same quantities are plotted in the panels (f)-(h) at 
the wave vector q = Qsdw. Here we choose the amplitude 
of the superconducting order parameter Ao = 2 meV. For 
the numerical purposes we also employ the broadening of the 
Green's function, 5 = 0.2 meV. 



In Fig. [ijc) and (d) wc present the imaginary and the 
real parts of Xo(Qafm,^) as a function of frequency uj 
at T = 1 K. In the non-SC state, the imaginary part 
is linear in lo at low frequencies which is a typical Lan- 
dau damping within the Fermi-liquid picture. In the SC 
phase, the imaginary part of the magnetic susceptibility 
becomes gapped. The magnitude of the gap, fig, is equal 
to 2Ao in the s-wave case. At larger frequencies Imxo in- 
creases slowly from zero. In comparison, for the d-wavc 





FIG. 2: (Color online) Calculated Fermi surface of the aig- 
band model. The position of the nodes of the dj.2_y2-wave 
[(a),(c)] and dxy-wstve [(b), (d)] superconducting gaps is de- 
noted by the solid (green) curves. The plus and the minus 
signs refer to the corresponding phases of the superconduct- 
ing order parameter. The states at the FS connected by the 
wave vectors Qafm or Qsdw are shown by the circles. 



try it is smaller than 2Ao. From Fig. ^a) and (b) one 
can notice that the Qafm wave vector connects parts 
of the FS where Ak = — Ak-t-q^FM but also some parts 
where Ak = +Ak+QAPM. For d2.2_j^2-wave supercon- 
ducting gap, there are four pairs of points of the first 
type and two pairs of points of the second type. Due to 
the smaller |Ak| + |Ak+QAFMl foi' the first process, as it is 
seen from Fig.[21Ia), the Imxo shows a discontinuous jump 
at rig . This is due to the change of sign in the anomalous 
coherence factor, C^^. The second process will give con- 
tribution at energies larger than fig due to larger value of 
I Ak| + 1 Ak+QAFM I there. Therefore, the net effect will re- 
sult in a discontinuous jump of Imxo at Q,g. Correspond- 
ingly, the real part will possess a logarithmic singularity 
as it is also seen in Fig. [Ijd). Within the RPA the for- 
mation of the pole (spin resonance) in the total magnetic 
susceptibility below fig is possible if Im[xo(q, ^)] = 
and simultaneously 1/U = Re[xo(q, w)]. Due to the log- 
arithmic character of the singularity this condition will 
be generally fulfilled for any small value of U which would 
give a position of the resonance exactly at or very close 
to fig. However, a small amount of impurities or disorder 
will smear the singularity out and suppress the resonance 
peak. In Naa;Co02-?/H20 the value of U should be rela- 
tively large to shift the position of the spin resonance 
towards energies smaller than fig and make it robust 
against impurity scattering. The calculated susceptibil- 
ity is shown in Fig.[lje) where we use Ures— 0-^79 eV. It 
is interesting to note that the resonance occurs for both 



d^2_y2- and 



_y2 +idxy-'vfave symmetries, however, the 



lA- 



k+Q 



at 



symmetries the lowest value of fig = jAj, 
the Fermi surface. Obviously for the non-s-wave symme- 



value of fig slightly differs. Note, for dxy-wave supercon- 
ducting gap the situation is opposite. From Fig. [2{b) 
one sees that in contrast to d2,2_j,2-wave case there are 
two pairs of points at the FS where Ak = — Ak+q^pM 
and four pairs of points where Ak = -f Ak+Q^FM • Here, 
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FIG. 3: Calculated temperature dependence of the Knight 
shift K{T) (a) and the spin-lattice relaxation rate 1/T\T 
(b) for the aig-band model. Note the logarithmic temper- 
ature scale in (b). Here, we assume the conventional BCS 
temperature dependence of superconducting gap, Ao(T) = 

Ao^yi-T/n. 



the fig is determined by the second process, thus there 
will be no logarithmic jump in Rexo a-t ^g- Of course it 
will occur at larger frequencies due to the first type of 
process but the resonance conditions will not be fulfilled. 
Therefore, we do not expect the spin resonance for the 
dajy-wave symmetry. 

The present value of Ures is of course too small to 
be the on-site Coulomb repulsion which is of the order 
of several electron volts. Therefore, the effective inter- 
action U entering our model ([T]) originates mainly from 
the Hund's exchange, Jh- In the lamellar sodium cobal- 
tate, the value of Jh is presently disputed and the lowest 
estimated value is of the order of leV— . It has been 
shown recently that even this value significantly affects 
the population of the aig and e'g orbitala^i^^i^. Taking 
this value into account, we assume U = aJn, where Jh 
is the mean-field value of the Hund's exchange and a is 
the coefficient that describes corrections beyond mean- 
field theory. One has to note that the larger value of U 
will lead to the SDW instability in our calculations. 

The situation changes for the wave vector Qsdw 
[Fig.[TJf)-(h)]. There is one striking difference in the low- 
energy behavior of Im[xo(QsDW, t^^)]- Namely, already 
in the normal state the scattering rate is non-linear for 
small Lu. It is obviously a consequence of the 2kF insta- 
bility and a resulting non-Landau damping at this wave 
vector. Furthermore, in the SC state the situation differs 
drastically with respect to Qafm- As one could see from 
Fig- and (d) there is equal number of contributions 
for which Ak = -Ak+QsD„ and Ak = -}-Ak+QsDw- As 
a result the discontinuity does not occur and the real part 
of xo is smaller in the superconducting state than in the 
normal state. Therefore, for reasonable values of U there 
is no resonance condition for xrpa [see Fig.[TJh)]. 

Generally, a formation of the resonance peak below Tc 
in the unconventional superconductors is a well-known 
consequence of the sign change of the superconducting 
order parameter. It has been originally discussed in rela- 
tion to the layered high- Tc cuprates^ and also recently 
has been used to explain the inelastic neutron scatter- 
ing results in heavy-fermion compound UPd2Al3^. In 
layered superconducting cobaltates the emergence of the 



resonance peak for several symmetries of the supercon- 
ducting order parameter has been analyzed within sim- 
ple single-band model^. In contrast to Ref. H^, we have 
found that the resonance peak (even within simple aig- 
band model) is very sensitive to the small variation of 
LZ-values and to disorder. As a result the resonance is 
confined to the wave vector Qafm and disappears for 

IQI < IQafmI- 

The temperature dependence of the Knight shift, 
K{T), and the spin-lattice relaxation rate, 1/TiT, is cal- 
culated according to the expressions: 

K{T) cx lim Rcx(q,c^-0), (4) 

1/TrT oc limiy^^^^^^. (5) 
q 

In Fig. Owe show both quantities as a function of temper- 
ature. In the normal state l/TiT increases with decreas- 
ing temperature that reflects the presence of the incom- 
mensurate antiferromagnetic fluctuations in this system. 
At the same time, the Knight shift is a constant which 
stresses that there are no small-q fluctuations. Below T^, 
both physical observables drop rapidly due to opening 
of the superconducting gap in the energy spectrum. As 
expected, the decrease is exponential for dx2_y2 + id^y- 
wave symmetry due to its nodeless character in E^. For 
(i2;2_j^2-wave symmetry the behavior of \/TiT and K{T) 
follows standard power-law temperature dependence due 
to the presence of the line nodes in the energy spectrum. 
In the next section we will compare our results to the ex- 
perimental data where we describe a more realistic model 
in application to the superconducting cobaltate. 

III. t29-BAND MODEL 

The oig-band model is, of course, oversimplified for 
describing the physics of Na^CoO 2 -2/112 since aig-e^ 
level splitting, (5e, is only 53 meV. As a result there is 
a substantial hybridization of the aig and the bands, 
completely neglected in the simple oig-band model. In 
particular, the bands may form hole pockets at the FS 
in addition to a large aig-pocket^. To take into account 
these details, we further analyze the magnetic response in 
the full t2g-band model including both aig and e'g cobalt 
states. 

The free electron Hamiltonian of the t2g-band model 
in a hole representation is given by 

where n^aa = dl^^d^aa, rfkacr {dlaa) ^^e annihilation 
(creation) operator for the t2g-iio\e with spin cr, orbital 
index a, and momentum k, is the hopping matrix 
element, e" is the single-electron energy, and ^ is the 
chemical potential. All of the in-plane hoppings and 
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FIG. 4: (Color online) Calculated results for the t2g-band 
model, (a) q-dependence of Im[xo(q, <^)1 and Im[xi?,p.4(q, <^)] 
at = 5 meV in the normal (non-SC) phase. The scattering 

wave vectors QaFM, QsDWI, QsDW2, Q'sDW; Qae', Qe', 

Qe" , and Qe'" are denoted by the arrows, (b) The calcu- 
lated Fermi surface with the corresponding scattering wave 
vectors. In (c) the band dispersion is shown where the bold 
(blue) curve denotes the topmost band used for the suscep- 
tibility calculations. A horizontal (green) line stands for the 
chemical potential. The panels (d)-(f) show imaginary and 
real parts of xo, and imaginary part of Xrpa at q = Qafm 
in the normal state and in SC state with various supercon- 
ducting order parameter symmetries. The imaginary parts 
of the bare and the total susceptibilities are plotted in the 
panels (g)-(h) and (i)-(j) at the wave vectors q = Qsdwi 
and q = Qsdw2 , respectively. Here we choose the ampli- 
tude of the superconducting order parameter Ao = 2 meV. 
For the numerical purposes we also employ the broadening of 
the Green's function, 5 = 0.2 meV. 



the single-electron energies were derived previously by 
us from the ab-initio LDA calculations using projection 
procedure and we use here the parameters for a;=0.33 
from Ref. Issl . To obtain the dispersion we diagonahze 
the Hamiltonian ([6]) calculating the chemical potential /i 
self-consistently. The resulting FS topology and energy 
dispersion arc shown in Fig. lUJb) and (c), respectively. 
The resulting dispersion and the FS replicate the corre- 
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FIG. 5; Calculated temperature dependence of the Knight 
shift K{T) (a) and the spin-lattice relaxation rate 1/TiT (b) 
for the t2g-ha,nd model. Note, in (b) the curve for [7=0.26 eV 
was scaled by a factor of 0.025. 



spending LDA ones^. 

Due to the non-zero inter-orbital hopping matrix ele- 
ments, aig and bands are hybridized. However, only 
one of the hybridized bands crosses the Fermi level thus 
making the largest contribution to the low-energy proper- 
ties of the system. We refer to this band as ek- Note, it is 
substantially different from the simple ai^-band. Later, 
this effective band ek will be used to calculate the dy- 
namical magnetic susceptibility with some effective on- 
site Coulomb interaction U. 

Present FS has more complicated structure in compar- 
ison to the aig-band model. First, states are present at 
the Fermi surface and strongly hybridize with aig states. 
At the same time, the "rounded hexagon" shape of the 
central part of the FS arises from the hoppings beyond 
nearest- neighbors included in the t2g-band model and 
neglected in aig-band model considered above. This re- 
sults in a number of additional scattering wave vectors 
as calculated from xo, see Fig. HJ^a). In particular, there 
are four scattering wave vectors connecting the e'g-e'g FS 
pockets [Qe', Qe", Qe'", and Qsdw2 = (0,0.495)], and 
also two scattering wave vectors connecting the aig-e^ 
FS pockets [Qae' and Q'sdwI- ^^'^ same time, these 
wave vectors also connect parts of the central aig FS 
pocket and the total magnetic susceptibility includes con- 
tribution from this scattering too. In addition, there 
are two wave vectors, [Qafm and Qsdwi = (0, 0.649)] 
which arise due to the curved form of the central aig 
FS pocket. The pronounced peaks at all these wave vec- 
tors are present in both the bare and the RPA magnetic 
susceptibility {U=0.15 cV). Again, similar to the aig- 
band model, the magnetic response is not dominated by 
the scattering at the commensurate wave vector Qafm- 
The overall picture of the magnetic response is consistent 
with the one presented in Ref. [l9l . 

In the non-SC phase and the SC phase with s-wave 
order parameter the behavior of x(q, w) at q = Qafm 
[see Fig. |l{d)-(f)] is similar to the one in the aig-band 
model. However, for the c?-wave symmetry of the order 
parameter, one finds that for uj > fig the states with 
equal signs of the superconducting order parameter (sec- 
ond type of the process ) contributes first, and the dis- 
continuous jump in Im[xo(QAFM, w)] occurs at higher 
energies. The particular form of the FS in the realis- 
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tic i2g-band model and more complicated band structure 
produce this effect. Therefore, the resonance peak in 
Im[xi?pA(QAFM, ^^)] may in principle still exist, how- 
ever, it occurs in a very narrow interval of the U values. 
This interval is determined by the resonance condition in 
the superconducting state and by the stability of a para- 
magnetic state above Tc. Here, we use ?7res=0.26 eV, 
which is more than twice smaller than in the oig-band 
model. 

Although the formation of the spin resonance is un- 
realistic for the antifcrromagnctic wave vector Qafm it 
may now occur at other wave vectors. In Fig. |4l^g)-(j) we 
present the imaginary parts of xoCQj'^) and XRPAi^l-,^) 
at QsDWi and at Qsdw2- Here, one notices the pro- 
nounced effects of the complicated i2g-band structure at 
high energies for the scattering at both wave vectors. 
Deviations from the linear-cj damping start already at 
low energies, smaller than Qg. For U = Ures the spin- 
resonance is present at Qsdwi for both c?-wave symme- 
tries. However, at Qsdw2 the resonance peak is present 
for (1^2 _y2 + idxy-'wavc symmetry only. Similar to the sit- 
uation with Qafm, this is due to smallness of the allowed 
U values. 

In Fig. [5] wc show the corresponding results for the 
1 /TiT and K{T) . Below superconducting transition tem- 
perature the behavior is very similar to the results ob- 
tained for the simple aig-band model. This is because be- 
low Tc the symmetry of the superconducting gap and its 
nodal structure determines the temperature dependen- 
cies of the l/TiT and the K{T) values. At the same time, 
notice the stronger AFM fluctuations in the normal state. 
For almost the same value of U as in Fig. [3] this is due 
to the larger density of states at the Fermi level (and the 
change of the Fermi velocity) than in the simple ai^-band 
model. Such a behavior is observed in the experimen- 
tal NQR datai^iiii^. It is interesting to note that with- 
out water the parent non-superconducting compound 
Nao.33Co02 shows much weaker AFM fluctuation&ii. In 
our theory the fluctuations occur for the parent com- 
pound too. It probably demonstrates a possible sig- 
nificance of the third dimension and, in particular, the 
bilayer splitting which may reduce the two-dimensional 
AFM fluctuations in Nao.33Co02. 

Note, the presence of the e'g pockets on the FS can also 
lift the degeneracy between the three c?-wave states. Since 
in the dx^-y^-wtcve SC state the e'g FS pockets arc fully 
gapped, the additional condensation energy is gained 
[compare the topology of the line-nodes in Fig. [2{a) and 
FS topology in Fig.m^b)]. For the c?a;,y-wave SC state this 
gain in energy will be smaller [compare Fig. [2jb) and FS 
in Fig.Hb)]. 

Presently, there is still a discussion on the details of the 
Fermi surface topology in the water intercalated cobal- 
tates. In particular, ARPES experiments do not ob- 
serve the Cg-pockets at the F S'^^'^^i'^^ . It has been shown 
that an inclusion of the electronic correlation within 
Gutzwillcr approximation may shift the e^-bands below 
the Fermi leve P^i^^ , although this conclusion has been 
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FIG. 6: (Color online) Calculated results for the t2g-band 
model with enlarged crystal field splitting, (a) q-dependence 
of the Im[xo(q, i^)] and the Im[xflPA(q, w)] at w = 5 meV 
in the normal (non-SC) phase. The scattering wave vec- 
tors Qafm, Qsdw, Q'sdwi and Qsdw are denoted by the 
arrows. The band dispersion is shown in the inset of (a), 
where the bold (blue) curve denotes the topmost band used 
for the susceptibility calculations, and the horizontal (green) 
line stands for the chemical potential, (b) The calculated 
Fermi surface with the corresponding scattering wave vec- 
tors. (c)-(e) The calculated imaginary (c) and real (d) parts 
of the Xo(Qafm, 'i'), and the imaginary part of Xrpa (e) in 
the normal and in the SC state with various superconducting 
order parameters. The same quantities are plotted in (f)-(h) 
at q = Qsdw. Here we choose the amplitude of the super- 
conducting order parameter Ao = 2 meV. For the numerical 
purposes we also employ the broadening of the Green's func- 
tion, 5 = 0.2 meV. 



challengecl22i^i2S. Another interpretation of this exper- 
imental result relays on the disorder introduced by Na. 
As it was shown within LDA, the scattering due to dis- 
order can destroy the small e^-pockets^. For the super- 
conducting polycrystalline samples, recent experiments 
indicate that the oxonium ions, HsO"'", may introduce 
additional dopanta^i^i^, or result in oxygen vacancies 
reducing Co oxidation stated. Though, this conclusion 
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T (K) T (K) 

FIG. 7: (Color online) Calculated temperature dependence 
of the Knight shift K{T) (a) and the spin-lattice relaxation 
rate 1/TiT (b) for the f2g-band model without e'g FS pockets. 
Note, in (b) the curve for [7=0.342 eV was scaled by a factor 
of 0.002. 

has been doubted by the NMR experiments^ which show 
the Co valence state is insensitive to hydration and de- 
pends on the Na content only. This was also confirmed 
later by the powder neutron diffraction^. 

In our study we further consider the t2g-band model 
with increased crystal filed splitting, Se—153 meV. This 
makes e'g band sink below the Fermi level, as it is seen 
in the inset of Fig. ^a.). The behavior of the dynamical 
spin susceptibility for U~Q.15 cV at w = 5 mcV pre- 
sented in Fig. EJa) shows more similarity to the simple 
oig-band model with additional features due to peculiar- 
ities ( "rounded hexagon" form) of the large FS pocket as 
shown in Fig. [6jb). The scattering is most pronounced 
at the wave vector Qsdw — (0,0.633). There is also 
intensive scattering at the wave vector Qsdw owing its 
appearance to the curved shape of the FS. 

Fig. [6l[c)-(e) and (f)"(h) displays the magnetic suscep- 
tibility at Qafm and at Qsdw, respectively. Contrary 
to both oig-band model and t2g-band model with e'g FS 
pockets, here wc observe a well-defined linear behavior of 
Im[xo(q7t^)] in the considered frequency range at these 
wave vectors. For the d-wave order parameter, the be- 
havior of the susceptibility resembles that in the t2g-band 
model with e'g FS pockets. Again one could find a nar- 
row range of parameters where the resonance peak exists, 
which we illustrate in Fig. [nje),(h) for [/,.es=0.342 eV. 

Similarly, the change of the FS topology does not in- 



fluence significantly the temperature dependence of the 
Knight shift and the spin-lattice relaxation rate above 
and below Tc- This is illustrated in Fig. [7] where we plot 
both quantities as a function of temperature. 



IV. CONCLUSION 

Our analysis of the dynamical spin susceptibility in 
application to the Na2;Co02-yH20 have shown that the 
magnetic response in the normal state is dominated by 
the incommensurate SDW fluctuations at large momenta 
close to Qafm- This is consistent with experimental 
NQR data which shows a pronounced AFM-like fluctu- 
ations in the temperature dependence of the spin-lattice 
relaxation rate. It is interesting to note that the pres- 
ence of the Cg-pockets at the Fermi surface is not affect- 
ing significantly this result. In the normal state we note 
the absence of ferromagnetic-like fluctuations. This ob- 
servation justifies our choice of spin-singlct order param- 
eter, because to induce the spin-triplet Cooper-pairing 
the fluctuations with small momenta are required. Be- 
low Tc our results for dx2^y2- or d^jy-wave (not shown) 
symmetries of the superconducting order parameter are 
consistent with experimental data which excludes node- 
Icss (1^2 _y2 4- idxy-wavc symmetry. We further stress that 
the resonance peak, predicted previously^ for the sim- 
ple flig-band model, is improbable for the realistic band 
structure of Na2;Co02-j/H20. Moreover, we find that even 
if present the resonance peak is confined to the AFM 
wave vector and disappears away from it. 



Acknowledgments 

We would like to thank A. Donkov, V. Yushankhai, A. 
Yaresko, and G.-q. Zheng for useful discussions, and I. 
Mazin for useful comments. M.M.K. acknowledge sup- 
port form INTAS (YS Grant 05-109-4891) and RFBR 
(Grants 06-02-16100, 06-02-90537-BNTS). I.E. acknowl- 
edge support form Volkswagen Foundation. 



* Electronic address: maxim@mpipks-drcsdcn . mpg. de' 

^ D. VoUhardt and P. Woelfle, The Superflmd Phases of He- 
lium 3 (Taylor and Francis, London, 1990). 

^ see H.F. Fong, P. Bourges, Y. Sidis, L.P. Regnault, J. 
Bossy, A. Ivanov, D.L. Milius, LA. Aksay, and B. Keimer, 
Phys. Rev. B 61, 14773 (2000) and references therein. 

^ K. Takada, H. Sakurai, E. Takayama-Muromachi, F. 
Izumi, R.A. Dilanian, and T. Sasaki, Nature 422, 53 
(2003). 

* G. Cao, C. Feng, Y. Xu, W. Lu, J. Shen, M. Fang, and 
Z.-a. Xu, J. Phys.: Condens. Matter 15, L519 (2003). 

^ B. Lorenz, J. Cmaidalka, R.L. Meng, and C.W. Chu, Phys- 

ica C 402, 106 (2004). 
« H.D. Yang, J.-Y. Lin, CP. Sun, Y.C. Rang, C.L. Huang, 



K. Takada, T. Sasaki, H. Sakurai, and E. Takayama- 
Muromachi, Phys. Rev. B. 71, 020504(R) (2005). 
R. Jin, B.C. Sales, S. Li, and D. Mandrus, Phys. Rev. B 
72, 060512(R) (2005). 

* A. Kanigel, A. Keren, L. Patlagan, K.B. Chashka, P. King, 
and A. Amato, Phys. Rev. Lett 92, 257007 (2004). 

^ T. Fujimoto, G.-q. Zheng, Y. Kitaoka, R.L. Meng, J. 
Cmaidalka, and C.W. Chu, Phys. Rev. Lett. 92, 047004 
(2004). 

K. Ishida, Y. Ihara, K. Kitagawa, H. Murakawa, Y. Maeno, 
C. Michioka, M. Kato, K. Yoshimura, K. Takada, T. 
Sasaki, H. Sakurai, and E. Takayama-Muromachi, J. Phys. 
Soc. Jpn. 72, 3041 (2003). 
^ G.-q. Zheng, K. Matano, R.L. Meng, J. Cmaidalka, and 



8 



C.W. Chu, J. Phys.: Condens. Matter 18, L63 (2006). 
Y. Ihara, H. Takeya, K. Ishida, H. Ikeda, C. Michioka, 
K. Yoshimura, K. Takada, T. Sasaki, H. Sakurai, and 
E. Takayama-Muromachi, J. Phys. Soc. Jpn. 75, 124714 
(2006). ' 

" C. Michioka, H. Ohta, Y. Itoh, and K. Yoshimura, J. Phys. 

Soc. Jpn. 75, 063701 (2006). 
^"^ M. Kato, C. Michioka, T. Waki, K. Yoshimura, K. Ishida, 

H. Sakurai, E. Takayama-Muromachi, K. Takada, T. 

Sasaki, Physica B 359-361, 482 (2005). 

Y. Ihara, K. Ishida, H. Takeya, C. Michioka, M. Kato, Y. 

Itoh, K. Yoshimura, K. Takada, T. Sasaki, H. Sakurai, and 

E. Takayama-Muromachi, J. Phys. Soc. Jpn. 75, 013708 

(2006) . 

" A. Tanaka and X. Hu, Phys. Rev. Lett. 91, 257006 (2003). 
" O.I. Motrunich and P.A. Lee, Phys. Rev. B 70, 024514 

(2004) . 

K. Kuroki, Y. Tanaka, and R. Arita, Phys. Rev. Lett. 93, 
077001 (2004). 

^® M.D. Johannes, I.I. Mazin, D.J. Singh, and D.A. Papacon- 
stantopoulos, Phys. Rev. Lett. 93, 097005 (2004). 

^° K. Kuroki, Y. Tanaka, and R. Arita, Phys. Rev. B 71, 
024506 (2005). 

M. Mochizuki, Y. Yanase, and M. Ogata, Phys. Rev. Lett. 

94, 147005 (2005). 
22 G.-q. Zheng, K. Matano, D.P. Chen, and C.T. Lin, Phys. 

Rev. B 73, 180503(R) (2006). 
2^ Y. Kobayashi, T. Moyoshi, H. Watanabe, M. Yokoi, and 

M. Sato, J. Phys. Soc. Jpn. 75, 074717 (2006). 
2* see I.I. Mazin and M.D. Johannes, Nature physics 1, 91 

(2005) and references therein. 

25 W. Higemoto, K. Ohishi, A. Koda, S.R. Saha, R. 
Kadono, K. Ishida, K. Takada, H. Sakurai, E. Takayama- 
Muromachi, and T. Sasaki, Phys. Rev. B 70, 134508 
(2004). 

W. Higemoto, K. Ohishi, A. Koda, R. Kadono, H. Saku- 
rai, K. Takada, E. Takayama-Muromachic, and T. Sasaki, 
Physica B 374, 274 (2006). 

2^ S. Florens and M. Vojta, Phys. Rev. B 71, 094516 (2005). 

2*^ H.W. Zandbergen, M. Foo, Q. Xu, V. Kumar, and R.J. 
Cava, Phys. Rev. B 70, 024101 (2004). 

2^ M. Mochizuki and M. Ogata, J. Phys. Soc. Jpn. 76, 013704 

(2007) . 

^° D.J. Singh, Phys. Rev. B 61, 13397 (2000). 

31 M.Z. Hasan, Y.-D. Chuang, D. Qian, Y.W. Li, Y. Kong, 
A. Kuprin, A.V. Fedorov, R. Kimmerling, E. Rotenberg, 
K. Rossnagel, Z. Hussain, H. Koh, N.S. Rogado, M.L. Foo, 
and R.J. Cava, Phys. Rev. Lett. 92, 246402 (2004). 



32 H.-B. Yang, S.-C. Wang, A.K.P. Sekharan, H. Matsui, S. 
Souma, T. Sato, T. Takahashi, T. Takeuchi, J.C. Cam- 
puzano, R. Jin, B.C. Sales, D. Mandrus, Z. Wang, and H. 
Ding, Phys. Rev. Lett. 92, 246403 (2004). 

33 T. Shimojima, K. Ishizaka, S. Tsuda, T. Kiss, T. Yokoya, 
A. Chainani, S. Shin, P. Badica, K. Yamada, and K. 
Togano, Phys. Rev. Lett. 97, 267003 (2006). 

3* S. Zhou, M. Gao, H. Ding, P.A. Lee, and Z. Wang, Phys. 

Rev. Lett. 94, 206401 (2005). 
35 M.M. Korshunov, I. Eremin, A. Shorikov, V.I. Anisimov, 

M. Renner, and W. Brenig, Phys. Rev. B 75, 094511 

(2007). 

3^ A.O. Shorikov, V.I. Anisimov, and M.M. Korshunov, cond- 

mat/0705.1408 (unpublished). 
3^^ H. Ishida, M.D. Johannes, and A. Liebsch, Phys. Rev. Lett. 

94, 196401 (2005). 
3** C. A. Perroni, H. Ishida, and A. Liebsch, Phys. Rev. B 75, 

045125 (2007). 

3^ A. Liebsch and H. Ishida, cond-mat/0705.3627 (unpub- 
lished) . 

^° see I. Eremin, D.K. Morr, A.V. Chubukov, K.H. Benne- 
mann, and M.R. Norman, Phys. Rev. Lett. 94, 147001 
(2005) and references therein. 

''I J. Chang, I. Eremin, P. Thalmeier, and P. Fulde, Phys. 
Rev. B 75, 024503 (2007). 

'*2 J.-X. Li and Z.D. Wang, Phys. Rev. B 70, 212512 (2004). 

''3 Y. Ihara, K. Ishida, K. Yoshimura, K. Takada, T. Sasaki, 

H. Sakurai, and E. Takayama-Muromachi, J. Phys. Soc. 
Jpn. 74, 2177 (2005). 

D.J. Singhl and D. Kasinathan, Phys. Rev. Lett. 97, 
016404 (2006). 

D.P. Chen, H.C. Chen, A. Maljuk, A. Kulakov, H. Zhang, 
P. Lemmens, and C.T. Lin, Phys. Rev. B 70, 024506 
(2004). 

C.J. Milne, D.N. Argyriou, A. Chemseddine, N. Aliouane, 
J. Veira, S. Landsgesell, and D. Alber, Phys. Rev. Lett. 
93, 247007 (2004). 
■'^ K. Takada, K. Fukuda, M. Osada, I. Nakai, F. Izumi, R.A. 
Dilanian, K. Kato, M. Takata, H. Sakurai, E. Takayama- 
Muromachi, and T. Sasaki, J. Mat. Chem. 14, 1448 (2004). 
P.W. Barnes, M. Avdeev, J.D. Jorgensen, D.G. Hinks, H. 
Glaus, and S. Short, Phys. Rev. B 72, 134515 (2005). 

I. R. Mukhamedshin, H. AUoul, G. CoUin, and N. Blan- 
chard, Phys. Rev. Lett. 94, 247602 (2005). 

5° L. Viciu, Q. Huang, and R. J. Cava, jCond-mat /0603600| 
(unpubhshed). 



